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Analysis of an object data set on a multi-dimensional space. 



The invention relates to a method of analyzing an object data set on a multi- 
dimensional space. 

Such an object data set represents one or more properties of the object to be 
examined. The object data set notably relates to the density distribution in the object to be 
5 examined; in that case the data values are the local density values of (a part of) the object to 
be examined. The data values may also relate, for example, to the distribution of the 
temperature or the magnetization in the object. The multi-dimensional space is usually the 
O three-dimensional space. The data values then relate to a volume distribution of the relevant 
|| property, for example, the density distribution in a volume of the object to be examined. The 
^ j 1 0 multi-dimensional space may also be two-dimensional. In that case the data values relate to a 
[II distribution of the relevant property in a plane through the object, for example the density 
h : distribution in a cross-section through the object. 

* n The object data set can be acquired in a variety of ways. The object data set 

m notably relates to a patient to be examined. Such an object data set can be acquired by means 
r is of various techniques such as 3D X-ray rotational angiography, computed tomography, 
^ magnetic resonance imaging or magnetic resonance angiography. 

if ^ 

The article "3D rotational angiography: Clinical value in endovascular 
treatment" in Medica Mundi 42 (1998) by J. Moret et al. concerns the determination of the 
section of a blood vessel. 

20 According to the known method a more or less tubular structure relating to the 

blood vessel of the patient to be examined is identified from the object data set. Such an 
identification is performed on the basis of differences between data values relating to the 
blood vessel and to the surrounding tissue. Such differences are realized by introducing a 
suitable contrast agent into the vascular system of the patient to be examined during the 

25 acquisition of the object data set. According to the known method a cutting plane is visually 
taken in the object data set so as to extend more or less transversely of the tubular structure. 
Subsequently, the diameter of the relevant blood vessel is calculated from the cross-section of 
the tubular structure along the cutting plane. 
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It is a drawback of the known method that it is difficult to orient the cutting 

plane so as to extend accurately perpendicularly to the tubular structure. This inaccuracy 

introduces serious errors in the calculated diameter of the relevant blood vessel. 

It is an object of the invention to provide a method wherein analysis of a 

tubular structure in the object data offers more accurate results than the known method. It is 

notably an object of the invention to provide a method whereby the diameter of such a 

tubular structure can be accurately determined. 

This object is achieved according to the invention by means of a method of 

analyzing an object data set in which a tubular structure occurs, wherein 

said object data set assigns data values to positions in a multi-dimensional 

space, which data values relate to an object to be examined, 

a starting position is chosen in or near the tubular structure, 

gradients to the surface of the tubular structure are derived at the area of the 

starting position, 

a cutting plane is derived through the starting position so as to have a direction 
as parallel as possible to the gradients in positions to the surface of the tubular structure and 
in a vicinity of the starting position, and 

a cross-section of the tubular structure is derived along the cutting plane. 

The object data set can be acquired, for example, by means of various 
techniques such as magnetic resonance angiography, computed tomography or 3D rotational 
X-ray angiography. Techniques of this kind produce an object data set with data values 
representing the structure of (a part of) the vascular system of the patient to be examined. 

The tubular structure concerns a part of the object data set in which the data 
values deviate from the surroundings of the tubular structure. The deviating data values lie, 
for example, in a selected range of data values which are larger than a predetermined 
threshold or smaller than a predetermined ceiling value. The positions in the multi- 
dimensional space for which the data values have such a deviating value are situated in an 
elongate, more or less tubular region. This tubular region may be shaped as a vacated or a 
filled tube. The tubular structure can be indicated in the object data set by means of 
segmentation techniques, thus identifying the positions with the deviating data values. 

The starting position indicates the position in which the cross-section of the 
tubular structure is desired. The starting position can be indicated, for example, by pointing it 
out in a rendition of the object data set on a display screen. 
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The gradients concern local differences between data values at the area of the 
surface. According to the invention the normal vector to the cutting plane is oriented exactly 
transversely of the gradients to the surface. The normal vector is oriented in such a manner 
that the normal vector is perpendicular to the gradients with a small tolerance only. The 
cutting plane can thus be oriented substantially parallel to the tubular structure. As a result, 
the diameter of the tubular structure can be accurately derived from the cross-section of the 
tubular structure. This diameter represents, for example, a useful technical result of the 
examination of the vascular system of the patient to be examined. It is notably possible to 
perform an accurate examination of stenoses or aneurysms. 



These and other aspects of the invention will be described in detail hereinafter 
on the basis of the following embodiments which are defined in the dependent Claims. 

The normal vector to the cutting plane is oriented, for example, by minimizing 
the weighted sum of the squares of the scalar products of the normal vector and the gradients 
to the surface of the tubular structure in a (small) vicinity of the starting position while 
varying the normal vector to the cutting plane. The weights in the weighted sum can be used, 
for example to make the effect of gradients in positions nearer to the starting point on the 
normal vector to the cutting plane greater than the effect of gradients in positions situated 
somewhat further from the starting position. 

Preferably, a local center of the tubular structure in the cutting plane is derived 
from positions in the cutting plane having data values in the selected range. This local center 
represents the center of the tubular structure in the cross-section along the cutting plane. For 
example, the local center is situated in a position which is situated at approximately the same 
distance from practically all positions in the cutting plane having data values in the selected 
range. It has been found that in many cases the vast majority of such positions having data 
values in the selected range and in the cutting plane represent locations at the edge of the 
tubular structure. Preferably, positions in which the gradients have a magnitude greater than a 
minimum value are used to derive the local center. Positions at the edge of the tubular 
structure are thus accurately selected. 

An even more accurate result for the local center is achieved iteratively by 
means of the method defined in Claim 4. Preferably, a subdivision into four equally large 
sectors is used. Positions outside the relevant tubular structure are excluded by determining 
the minimum distance of positions having data values in the selected range in each sector. 
When different minimum distances from the current center position of the local center occur 
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in different sectors, the current center position will not be situated exactly at the center of the 
tubular structure, A more accurate estimate of the center position is obtained by shifting the 
center position slightly in the direction of the sector in which the largest minimum distance 
occurs. Preferably, the center position is shifted over a distance corresponding to half the 
5 difference between the minimum distances in oppositely situated sectors. 



the tubular structure. To this end there is derived a subsequent position which is shifted in the 
direction of the normal to the cutting plane in the starting position. A next cutting plane 
through the subsequent position is then derived in conformity with Claim 1 . In the next 
10 cutting plane a local center of the cross-section through the tubular structure is determined 
again. Successive local centers are derived along the tubular structure by repetition. The 
magnitude of the shift along the tubular structure can be chosen by the user. The tracking of 
O the tubular structure will be more exact as the shifts for deriving each time a next subsequent 
m position are smaller. Preferably, a series of diameters is calculated for the series of cross- 
J"! 15 sections through the tubular structure thus obtained along respective cutting planes. This 
f II offers a thorough insight into the variation of the diameter of the tubular structure along the 
JH longitudinal axis of the structure. This result constitutes a useful aid for examining the extent 
* of stenosis in the vascular system of the patient to be examined. 



a diameter which is significantly larger than that of the inlet and outlet ducts. Such a situation 
occurs, for example, in the case of an aneurysm. The actual aneurysm then corresponds to the 
reservoir whereas the inlet and outlet ducts correspond to the blood vessels feeding and 
draining the blood to and from the actual aneurysm. The method disclosed in Claim 6 is 

25 preferably used in such a situation. According to the invention the inlet and outlet points can 
be determined by tracking the inlet and outlet ducts in conformity with Claim 5. Preferably, 
such tracking of the inlet and outlet ducts is observed on a monitor by the user who can 
terminate the tracking operation as soon as the connections to the reservoir are reached. 
Subsequently, a central position is determined in the reservoir. For example, a part of the 

30 object data set which contains the reservoir and a small (relative to the reservoir) number of 
positions outside the reservoir is selected. It has been found that the center of gravity of the 
data values in the selected part of the object data set suitably corresponds to the center of the 
reservoir. The cutting plane through the inlet and outlet points and through the center of the 
reservoir separates the reservoir from the inlet and outlet ducts. The volume of the reservoir 



According to the invention it is also possible to track the longitudinal axis of 



r:2o 



The invention is preferably used for tracking a tubular structure which 
includes a reservoir whereto inlet and outlet ducts are connected. The reservoir typically has 
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to the side of the cutting surface which faces the central position accurately represents the 
effective volume of the reservoir. When the reservoir with the inlet and outlet ducts 
represents an aneurysm and the connected blood vessels, the effective volume of the reservoir 
constitutes an accurate result concerning the size of the aneurysm; this is of importance for 
5 the further treatment of the aneurysm. 

Preferably, the method is carried out by means of a suitably programmed 
workstation. For example, a computer program with instructions for carrying out the steps as 
defined in the Claims 1 to 7 is loaded into such a workstation. 

Claim 7 defines a further application of the invention for a tubular structure 
1 0 with a reservoir whereto inlet and outlet ducts are connected. The inlet and outlet points are 
determined again by tracking the inlet and outlet ducts in conformity with Claim 5. 
According to this application of the invention a connection duct is interpolated between the 
inlet point and the outlet point and through the reservoir. Subsequently, the connection duct is 
isolated from the reservoir. It has been found that the remainder constitutes a good 
15 approximation of the effective reservoir which can be suitably used for determining the size 
of the aneurysm. 

These and other aspects of the invention will be elucidated, by way of non- 
limitative example, with reference to the following embodiments and the accompanying 
drawing; therein: 

20 Fig. 1 shows an example of an object data set with a tubular structure in which 

the application of the invention is indicated, and 

Fig. 2 shows an example of a flow chart of the method according to the 

invention. 

Fig. 1 shows an example of an object data set containing a tubular structure in 
25 which the application of the invention is indicated. Fig. 2 shows an example of a flow chart 
of the method according to the invention. Fig. 1 notably shows the tubular structure, in this 
case a rendition of a part of the vascular system of the patient to be examined, which has 
already been segmented from the possibly larger object data set. This is represented by the 
step 10 in Fig. 2. The segmentation consists, for example, in that all data values are set to a 
30 fixed value in positions in which the original object data set contains data values outside the 
selected range. The vascular system includes an aneurysm 3 which acts as a reservoir, a 
blood vessel 2 which transports blood to the aneurysm 3 and a blood vessel 4 via which 
blood is discharged from the aneurysm 3. The blood vessels 2 and 4 thus act as an inlet duct 
and an outlet duct. The starting point (BP) is selected by the user, for example by pointing it 
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out on the display screen of a workstation by means of a mouse and/or a keyboard and a 
cursor. Around the starting point (BP) there is taken a sphere (Bl) as indicated in step 20. 
The gradients (gy) of the data values are calculated at points on the edge of the tubular 
structure in the sphere (Bl); this is indicated in step 30. More explicitly the gradients are 
calculated as follows: 



9u = 



dXj 
DP 

dD 

K dZ U 



where (xj,yj,Zj) is the position at issue at the surface of the tubular structure and 



D indicates the data-values of the data-set which represent for example the pixel-values or the 
voxel- values of the object to be examined. Consequently, the vector gy points transverse to 
the surface of the tubular structure. Subsequently, in step 40 the normal vector h x is 

10 calculated by means of a minimization process, so that the sum J] Wj (n } • g Xj ) 2 is 

minimum. The weight factors wj are to be adjusted by the user; preferably, the weight factors 
decrease as the distance between the position in which the gradient gij is calculated and the 
starting point BP is larger. The cutting plane SN1 has a normal vector and extends through 
the starting point BP. The cross-section through the tubular structure 1 along the cutting 

15 plane SN1 is taken in step 50. In step 51 the local diameter d of the tubular structure 1, i.e. of 
the blood vessel, is derived from this cross-section. Furthermore, in step 52 the local center 
position of the tubular structure is derived from the cross-section by executing the iterative 
procedure in conformity with. Claim 4. 

In order to track the tubular structure as from the starting point BP, in step 60 

20 the subsequent position VP is derived in the direction of the normal vector h x which is locally 
oriented accurately along the axis of the blood vessel. The magnitude of the shift is adjusted 
by the user in dependence on the degree of meandering of the tubular structure. Fig. 1 shows 
a subsequent position VP which is reached after a large number of shifts along the axis of the 
blood vessel 2. When the subsequent position is not yet the end position (EP) desired by the 

25 user, the steps 20, 30, 40, 50, 51, 52 and 60 are repeated. Thus, the axis along the blood 
vessel and the diameter along the axis are thus determined. 

For example, the end point EP is the inlet point TP from the blood vessel 2 to 
the aneurysm 3. The outlet point AP is reached by tracking the blood vessel from a different 
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starting point BP f on the blood vessel 4 by means of the procedure shown in Fig. 2. For 
example, the user himself can terminate the tracking of the blood vessel when, on the basis of 
the anatomical insight of the user, the inlet point or outlet point has been reached. It is also 
possible to terminate the tracking of the blood vessel when the local diameter of the blood 
5 vessel suddenly increases strongly. A connecting line 1 is drawn through the inlet point and 
the outlet point (TP, AP). Furthermore, the central position Z of the aneurysm 3 is calculated 
as the center of gravity on the basis of the data values in the region G. The region G is, for 
example, a sphere around the aneurysm 3. It has been found that suitable results are achieved 
for the center position Z when it is ensured that the region G does not include too many 
10 points in the segmented data set which lie outside the aneurysm. The normal vector m is 

directed opposite to the perpendicular from Z to the connecting line 1. Subsequently, a cutting 
plane SN3 is determined through the connecting line and with the normal vector m . It 
appears that this cutting plane SN3 accurately separates the inlet and outlet ducts, in this case 
being the blood vessels 2 and 4, from the aneurysm 3. 
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1 . A method of analyzing an object data set in which a tubular structure occurs, 

wherein 

said object data set assigns data values to positions in a multi-dimensional 
space, which data values relate to an object to be examined, 
5 a starting position is chosen in or near the tubular structure, 

gradients to the surface of the tubular structure are derived at the area of the 
starting position, 

a cutting plane is derived through the starting position so as to have a direction 
in ®s parallel as possible to the gradients in positions to the surface of the tubular structure and 

10 in a vicinity of the starting position, and 
\! - a cross-section of the tubular structure is derived along the cutting plane. 

Vr 2. A method as claimed in Claim 1 , wherein the direction of the cutting plane is 

O derived by minimization of the sum of the squares of the scalar products of the normal vector 
1 5 to the cutting plane and the gradients. 

I- 3. A method as claimed in Claim 1, wherein a local center of the tubular structure 

in the cutting plane is derived from positions in the cutting plane having data values within a 
predetermined range, such a range concerning notably data values higher than a threshold 
20 value or data values smaller than a ceiling value. 

4. A method as claimed in Claim 3, wherein the cutting plane is subdivided into a 

plurality of sectors and 

a center position of the local center of the cross-section of the tubular structure 
25 in the cutting plane is estimated, 

for the respective sectors a minimum distance is derived between positions 
having data values in the preselected range and a current estimated center position, 

a maximum of the minimum distances in the respective sectors is determined, 
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a new estimate of the center position is derived, the new estimate of the center 
position being shifted, relative to the current estimate, in the cutting plane in the direction of 
the sector in which the maximum of the minimum distances is situated and over a distance 
which is dependent on the differences between positions in which the minimum distances 
occur in oppositely situated sectors, 

the steps b, c and d are repeated, if necessary, the new estimate of the center 
position being used each time as the current estimate. 



5. A method as claimed in Claim 3, wherein 

a subsequent position which has been shifted, relative to the starting position, 
in the direction of the normal to the cutting plane through the starting position is determined, 
and 

a subsequent cutting plane through the subsequent position is derived, the 
direction of said subsequent cutting plane being as parallel as possible to the gradients in said 
subsequent cutting plane and to the surface of the tubular structure. 



6. A method of analyzing an object data set in which a tubular structure occurs, 

wherein 

said object data set adds data values to positions in a multi-dimensional space, 

an inlet duct, an outlet duct and a reservoir occur in the tubular structure, the 
inlet and outlet ducts being connected to the reservoir, and 

an inlet point is determined in the inlet duct at the area of connection of the 
inlet duct to the reservoir, 

an outlet point is determined in the outlet duct at the area of connection of the 
outlet duct to the reservoir, 

a central position of the reservoir is determined, and 

a cutting plane through the inlet and outlet points is derived , the normal to the 
cutting plane extending along the perpendicular from the central position to the line through 
the inlet and outlet points. 



7. 

wherein 



A method of analyzing an object data set in which a tubular structure occurs, 
said object data set adds data values to positions in a multi-dimensional space, 
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an inlet duct, an outlet duct and a reservoir occur in the tubular structure, the 
inlet and outlet ducts being connected to the reservoir, and 

an inlet point is determined in the inlet duct at the area of connection of the 
inlet duct to the reservoir, 
5 an outlet point is determined in the outlet duct at the area of connection of the 

outlet duct to the reservoir, 

a connecting duct is interpolated between the inlet and outlet points and 
through the reservoir, and 

a remainder is derived as the difference between the reservoir and the 
1 0 interpolated connecting duct. 

8. A computer programme for analysing an object data set in which a tubular 

structure occurs, 

- the object data set assigning data values to positions in a multi-dimensional space, which 
1 5 data values relate to an object to be examined, 

— the computer programme comprising instructions for 

- choosing starting position in or near the tubular structure 

- deriving gradients to the surface of the tubular structure are derived at the area of the 
starting position 

20 - deriving a cutting plane through the starting position so as to have a direction as 

parallel as possible to the gradients to the surface of the tubular structure in positions 
in the vicinity of the starting position and 

- deriving a cross section of the tubular structure along the cutting plane. 
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ABSTRACT: 



The method relates to analysis of a tubular structure 1 such as a rendition of 
the vascular system of a patient. Gradients gy to the surface of the tubular structure 1 are 
computed in a starting point (BP). Subsequently, a normal h which extends essentially 
perpendicularly to the gradients is derived and a cross-section is taken along a cutting plane 
having h as its normal. A subsequent position (VP) is taken by performing a small shift in 
the direction of n . The axis of the vessel is tracked by repeating the procedure. 

Fig. 1 
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America filed by me on the same subject matter having a filing date before that of the application(s) of which priority is 
claimed: 



PRIOR FOREIGN/PCT APPLICATION(S) AND ANY PRIORITY CLAIMS UNDER 35 U.S.C. 119: 



COUNTRY 


APPLICATION NUMBER 


DATE OF FILING 
DAY, MONTH, YEAR 


PRIORITY 
CLAIMED UNDER 
35 USC 119 


Europe 


99203207.8 


1 October 1999 


YES 
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Combined Declaration For Patent Application and Power of Attorney (Continued) 
(includes Reference to PCT International Applications) 


Attorneys Docket Number 

PHN 17.678 US 


POWER OF ATTORNEY: As a named inventor, 1 hereby appoint the following attorney(s) abnd/or agent(s) to prosecute this 
application and transact all business in the Patent and Trademark Office connected therewith. (List name and registration number) 


Jack E. Haken, Rea. No. 26,902 
Michael E. Marion, Reg. 32,266 
Edward M. Blocker, Reg7T5o~"30^45 


Direct Telephone Calls to: 
(name and telephone number) 
(914)332-0222 



CO 



FULL NAME 
OF 

INVENTOR 



FAMILY NAME 
BRUIJNS 



FIRST GIVEN NAME 
Johannes 



SECONDE GIVEN NAME 



20 
1 



RESIDENCE 
& 

CITIZENSHIP 



CITY 

Eindhoven 



STATE OR FOREIGN COUNTRY 
The Netherlands j ^ 



COUNTRY OF CITIZENSHIP 
The Netherlands 



POST OFFICE 
ADDRESS 



POST OFFICE ADDRESS 
Prof. Holstlaan 6 



CITY 

5656 AA Eindhoven 



STATE & ZIP CODE/COUNTRY 
The Netherlands 



hereby declare that all statements made herein of my own knowledge are true and that all statements made on information and belief 
are believed to be true: and further that these statements were made with the knowledge that willful false statements and the like so 
made are punishable by fine or imprisonment, or both, under section 1001 if Title 18 of the United states Code, and that such willful 
false statements may jeopardize the validity of the application or any patent issuing thereon. 




SIGNATURE OF INVENTOR 201 
D4TE 26 April 2001 ~ 



U.S. DEPARTMENT OF COMMERCE- Patent and Trademarks Office 
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